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ABSTRACT 

The short time evolution of three dimensional small perturbations is studied. Exhibit- 
ing spectral asymptotic stability, thin discs are nonetheless shown to host intensive 
hydrodynamical activity in the shape of non modal growth of initial small perturba- 
tions. Two mechanisms that lead to such behavior are identified and studied, namely, 
non-resonant excitation of vertically confined sound waves by stable planar inertia- 
coriolis modes that results in linear growth with time, as well as resonant coupling 
of those two modes that leads to a quadratic growth of the initial perturbations. It 
is further speculated that the non modal growth can give rise to secondary strato- 
rotational instabilities and thus lead to a new route to turbulence generation in thin 
discs. 
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1 INTRODUCTION 

Reverting to Rayleigh's stability criterion Keplerian accretion discs have been deemed centrifugally stable. This has lead to 
focusing the attention on the magneto-ro t ationa l instability (MRI) which is excited due to the freezing of magnetic field 
lines into fluid elements iBalbus fc Hawlevl (|199ll ). Indeed, the MRI has been accepted since as a main turbulence generator 
and has consequently been widely considered to hold the key to the solution of the angular transfer problem in accretion 
discs. However, the MRI as a driver of turbulence and angular momentum tr ansport has recently been s hown to suffer from 
non-trivial difficulties in the sheer use of the shearing box in the simulat ions Regev fc Umurhan ( 200S ). Bodo et al. (2008)1 
as wel l as in doubts re g ardin g the numerica l conv ergence and resolution Lesur fc Longarettil (|2007l ). Fromang fc Papaloizou 
1 2007T ) . iFromang et al.l (|2007l ). iPessah et al.l (|2007l )]. In addition, it has recently been shown that if the thin disk geometry is 
taken into account both the gro wth rates as well a s the n umber of unstable MRI m odes are greatly reduced and are decreasing 
functions of the disk thickness [Coppi fc Kevesl (|2003h . iLiverts fc Mondl l|2009h ]. In parallel, various works have indicated 



over the last decade that pure hydrodynamical processes are not as irrelevant as possible sources for enhanced transport 
coefficients in thin accretion discs as has been hastily assumed before. In particular, two topics have caught the attention of 
researchers. In the first one the destabilizing effect of a stable stratification on a centrifugally stable rotating fluid has been 
exami ned. The result i ng ins t ability, termed the st r ato-rotational instabili t y (SR I), has been investigated for Taylor-Couette 



flows 



Boubnov et~aD (|l995h . lYavneh et aD |200lh . Ishalvbkov fc Rwdiger 



Dubrulle et al 



(2005)] or for Taylor-Cou e tte-like flows wi t hin th e 



( 2005t h lUmurhanl feood h iTevzadze et al.l HqO 



shearing sheet approximation that models rotating discs 
Common to those works is: 1. Rigid-like boundary conditions at two radial locations, 2. Radial pressure distribution plays 
a crucial role in the onset of the SRI under those conditions, 3. Ignoring the vertical boundaries of the fluid domain, and 
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4. Incompressibility within the Boussinesq approximation. We will return later on to these points (actually 1 and 2 are 
two sides of the same phenomenon) but let us first move on to the second route that may lead to pure hydrodynamical 
turbulence. The latter is motivated by succe s sful attemp t s by fl uid dynamics to expla i n the transition to turbulence in 



spectrally stable flows Boberg fc Brosa ( 19881 ). Grossmann ( 2000l ) . Schmid fc Henningson (f2000hl . Thus, re n dering the linear 



operator non normal, the flow shea r may give rise to trans i ently growing, perturba tions [Gustavssonl ( 199ll ) , iButler fc Farre' 



( 1992h . [Reddv fc Henningson! l| 19931 ). iTrefethen et all (|l993l ). lcriminale et all l| 19971 ). and the recent review by 
whose energy is eventually redistributed among different length-scales due to nonlinear effects. Indeed, recent calculations in 
real thin disc geometry have demonstrated the efficiency of such processes to significantly amplify initially sma ll perturbations, 
and th u s to generate in t ensive hydrodynamical activity in the otherwise centrifugally stable Keplerian discs 



Umurhan et al 



( 200d ) . iRebusco etafl (|2009l )]. Such non modal growth o f pertu rbations has also been obtained for ballooning modes in 



magnetized rotating plasmas in tokamaks [Chun fc Hameiril II 199(3 )]. The two routes described above are shown in the current 



work to diverge from a single unified comprehensive stability analysis in real thin disc geometry. It is thus shown first that if the 
proper physical conditions as well as the true thin disc geometry are taken into account, no spectrally unstable strato-rotation 
modes exist. Indeed, the lack of radial rigid boundaries, and the negligible role of the radial pressure gradients in establishing 
the steady-state rotation stand in stark distinction from the conditions that give rise to the SRI in Couette flows. Furthermore, 
the small thickness of the disc combined with its steady-state supersonic rotation results in vertical sound crossing time that 
is of the order of a rotation period. Consequently, compressibility effects cannot be ignored so that vertically propagating 
sound waves play a major role in the dynamical response of thin discs to small perturbations. 

Thus, (temporary) disappearing from the thin discs global scene (they will resurface later on), the SRI leaves the stage 
entirely to algebraic instabilities as a primary source of hydrodynamical activity. As will be shown there are two mechanisms 
that give rise to such non- exponential dynamical response. One is the non-resonant coupling between the two modes of wave 
propagation in the system (the inertia-coriolis waves and the sound waves). Such coupling is a direct result of the rotation 
shear that renders the linear operator non normal. Due to such interaction, sound waves are non-resonantly driven by the 
inertia-coriolis modes and their amplitude grows linearly in time. The second source of non-exponential perturbation growth 
is the resonant interaction between the above mentioned modes. In such scenario, sound waves are resonantly driven by the 
inertia-coriolis modes and due to the shear their amplitude grows quadratically in time. 

Finally, the two routes converge again by the tendency of the algebraically growing perturbations to develop sharp radial 
pressure gradients due to the rotation shear. As such gradients are essential ingredients in the onset of the SRI, it is conjectured 
that the latter may be excited as a secondary instability on top of the non modal growth. 

The paper is organized as follows. The dimensionless governing equations and their approximation to leading order in 
small aspect ratio of the steady-state disc are presented in the next Section. Section 3 describes the principal hydrodynamic 
modes of waves which can propagate in Keplerian discs. Non-resonantly as well as resonantly driven sound waves by inertial 
waves are considered in Sections 4 and 5, respectively. Conclusions are presented in section 6. Appendix A contains the 
principal boundary-value problems and methods of their solutions. 



2 THE GOVERNING EQUATIONS 



The dynamical response of radially and axially stratified thin rotating fluid discs to small perturbations is considered. For 
that purpose it is c onvenient to transform all the physical variables to non-dimensional quantities by using the following 
characteristic values [Shtemler et al.l (|2007l ) . (2009)]: 



v. = n.r., u = -, 



<E>, = V* 



P, = K(m»n* 



(1) 



Here = (GMc/rJ) 1 ' 2 is the Keplerian angular velocity of the fluid at the characteristic radius r, that belongs to the 
Keplerian portion of the disc ; G is the gravitational constant; M c is the total mass of the central object; is the characteristic 
value of the gravitational potential; m* and n, are the characteristic mass and number density; rm and m are the ion mass 
and number density; K is the dimensional constant in the polytropic law P = KrC , 7 is the polytrophic coefficient, (7 = 5/3 
in the adiabatic case). The resulting dimensionless dynamical equations, expressed in standard cylindrical coordinates {r, <fi, z} 
are: 

(2) 
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dt dr r dd> dz 
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1 1 dP 
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(3) 
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dn 1 drnVr 1 driVi dnV z , . 

^7 H 5 1 xT~ + — 5 — = °' ^ 5 - ) 

at r or r o(p oz 

P = n 1 . (6) 

Here the Mach number is defined as Ms = V,/Cs*, the characteristic sound velocity is Cs* = \J P*/(m»n*, t is time, 
V^>, 14} is the plasma velocity, P and n are the pressure and density current, and $ is the gravitational potential due to 
the central object. Note that a preferred direction is tacitly defined here, namely, the positive direction of the z axis is chosen 
according to positive Keplerian rotation. 



A common property of thin Keplerian discs is their highly compressible motion with large Mach numbers [ iFrank. King, fc Raine 

(2Q03)]: 



<1, 



(<) 



1 _IU 

Ms _ 6 ~~ ~ 

where the characteristic effective semi-thickness H„ of the equilibrium disc is such that the disc aspect ratio e equals the 
inverse Mach number (for the magnetic-field- free systems under consideration, the valueff* should be fixed and known). The 
smallness of e means that dimensionless axial coordinate is also small, i.e. z/r* ~ e (| z |^ H t ), and consequen tly it is 



conven i ent to introdu c e the following rescaled v alues in order to further a pply the asymp t otic e xpansions in e [see iRegev 
' 19831 ) . iQgilviel l| 19971 ) . iKluzniak fc Kital (|2000l ) . IRegev fc Gitelmanl (|200d ). lshtemler et alJ |2007h . (2009)]: 



c = - 

e 



H ±1 = ±Zl 
e 



(8) 



where the vertical velocity is scaled along with the disc height in e as follows from the principle of the least degeneracy of the 
problem for small e. 

The set of Eqs. ©-© is complemented with the conventional dynamic and kinematical boundary conditions on the 
interfaces, according to which the disc edges z = H±(t,r, <f>, z) are determined self-consistently. However, to simplify the 
problem, it is assumed that the disc edges are given by z — ±H* = ±const, where H t the characteristic effective semi-thickness 
of the equilibrium disc. Consequently, the kinematic boundary conditions are not needed, while the dynamic conditions 
P — n 1 — at the interfaces are modelled by: 

n = for C = ±1- (9) 

Using of the model conditions instead of the true boundary conditions at the interfaces indeed leads to the loss of the 
corresponding interface equations. However, this does not lead to a change in the essential nature of the dynamical behavior 
of the disc, and the kinematic boundary conditions provide separate equations for calculating the perturbed of the top and 
bottom interfaces of the disc. 

The gravitational potential in Eqs. ([2J and ((4]) may be expanded now in terms of the stretched axial variable to yield: 



--. ^ 2) (r,0 = ^ 1 , r>l»e. 
r 2r 6 



$(r, C) = $ (0) (r) + e 2 $ (2) (r, C) + 0(e 4 ) , $ <0) (r) 

Representing the solution of Eqs. ([I])-© as a sum of the equilibrium state and its perturbations 
P = P + P', n = n + n, V r = V T , V^ = V <t> + V^, Vk = Vl, 
the steady-state solution is given by 

P(r,0 = S 7 , ^'0=( 1 ^ i )^ T ( 1 -r^) 7iT , ^(r)=rfl(r), O(r) = r^ 2 . 

Substituting Eqs. (fT))- (|12[) into (2j|- ((9jl , linearizing about the steady-state solution, and employing the relation P' = c|n' 
yield the following linear set of equations for the perturbations to leading order in e: 



(10) 
(11) 
(12) 



dv' a( .dv; 



mv; = o, 



dnV' z 
dt 



n(r 



*n 2 (r)-£n' + n 2 (r)c 2 s — = 0, 



On' , dn' dnVt 

9T + " W W + ^ 



d<p 7-1 lj rfC 

1 drnVr _ 1 dnV£ 
r dr r d(f> 



where 

X =2 , c 5 (C) 

r dr 



Cs(r,0 _ 7-1 
ST2 2 (r) 2 



(i-C 2 ), 



(13) 
(14) 
(15) 
(16) 

(17) 
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C s = {dP/dnf/ 2 is the dimensionless equilibrium sound velocity; \ i s the epicyclic frequency. For Keplerian discs the 
epicyclical frequency is equals to Keplerian frequency \ — f2(r). 

Following Eq. (7), Eqs. (13)-(16) are supplemented with the following boundary conditions: 

n = for C = ±1- (18) 



3 THE MODES OF WAVE PROPAGATION 



3.1 Planar inertia-Coriolis waves 

We start the solution of the linearized system by noticing that Eqs. (|13|l and (|14p for V r ' and are decoupled from the rest 
of the system due to the thin disc approximation. Furthermore, the sub-set (|13J1 - (|14[ 1 represents a set of ordinary differential 
equations in time and only in the spatial variable 4>, while both the axial as well as the radial coordinates are parameters. It 
should be stressed that this is not an extra approximation but rather follows from the special geometry under investigation. 
In particular, it is noticed that that sub-set does not contain any pressure effects. The absence of pressure gradients in the 
horizontal dynamics hints also to the absence of the SRI in thin discs. Thus, assuming the following form for the solution of 
the sub-set (fl3 ]> - (fT3 ]> : 

where m — 0, ±1, ±2, ... is the azimuthal wave number, ui is the frequency of perturbations. The reason for transforming to 
the scaled mass flux variables defined in Eqs. (|19p is that, as will be shown below, most of the resulting equations do not 
depend on the radial coordinate. The corresponding dispersion relation is given by: 



{v' r , v'$y = exp(iujt — im(f>){vr, + c.c.. 



(19) 



to = m£l(r) ± x( r )- 



(20) 



This represents, not surprisingly, the stable epicyclic oscillations in the disc plane. As r and C, are mere parameters each ring 
r = const, ( = const vibrates independently in its own plane. If viscous stresses are taken into account, an axial profile is 
imposed due to the mutual shear stress es between the r i ngs a nd e ach entire cylindrical shell vibrates independently. This is 
indeed the case that has been solved in [Umurhan et all (fcood ) and lRebusco et all |2009h ]. Below x = ^( r ) f° r the Keplerian 
rotation will be used for simplicity. 

As the axial and radial coordinates play the role of passive parameters, the eigenfunctions of the inertia-Coriolis modes 
are determined up to arbitrary amplitude A(r, £): 

.1 



As an example the following form of the amplitudes is considered 
A(r,() = F(r)G(C), 



(21) 



(22) 



where the F(r) and G(£) are arbitrary functions. Any special form of functions A(r, £), or F(r) and G(£) represent some 
given set of initial conditions. That special form of the amplitude may be utilized to satisfy various regularity conditions of 
the solutions near the disc's edges. The regularity problem will be addressed in more details in sections 4 and 5. 

Finally, it is of interest to follow the evolution in time of the radial dependence of the perturbations. Thus, defining a 
generalized radial wave number by, say, k r = idlnVr/dr, it is readily seen that it depends on time as k r — — tXdfl/dr with 
the constant A. Thus, even initial radially uniform per turba tion will eventually deve lop strong radial dependence due to the 



rotation shear [see also lGoldreich fc Lvnden-BelJ l| 19651 ). and lBalbus fc Hawlevl (| 199 1 )] . As a result, in a general case an initial 



value problem has to be solved rather than a spectral one. In the current supersonically rotating thin disc case, however, due 
to the absence of radial partial derivative from the linearized sub-system (which is due to the insignificant role played by the 
radial steady-state as well as perturbed pressure gradients) a spectral problem may be solved in which the generalized time 
dependent radial wave number may be inferred from the solution. 



3.2 Vertical sound waves 

For frequencies that are different from w = mQ(r) ± x( r ) m l|20[l both perturbed components of the horizontal mass flux are 
zero, i.e. v r = and v^, = 0. As a result, Eqs. (|15|) and (I16[) describe the perturbations in the axial velocity and the number (or 
mass) density. Since v r = the radial derivative is once again dropped from the problem (from Eq. (|16|) V Thus, the following 
normal-wave solution is assumed: 

{v' z , n'} = exp(iu;i — im<j>){v z , h} + cc, (v' z = -pp) ■ (23) 
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Writing w = ASl(r) where A is a constant (an assumption that will be justified a posteriori) and substituting (|21l) - (|23[) into 
JTSJ -JTBJ yields: 

i(A-m)«.-^2^A + 4(0^=0, (24) 

z(A-m)n+%=0. (25) 
Finally, Eqs. (|24[) and (|25|) result in the following second order ordinary differential equation: 

In general, the solutions of Eq. (|26|l . complemented by the boundary conditions h — at the disc edges £ = ±1, may be 
expressed in terms of hypergeometric functions. However, for 7 = 5/3 the solutions become remarkably simple and may be 
expressed in terms of the following symmetric and anti-symmetric eigenfunctions (Appendix A): 

a k = 2k + l, n = N 2k +i cos(2fc + 1)0, k = 0, ±1, ±2, . . . , (27) 

a fe = 2k, n = N 2k sm2k6, k = ±1, ±2, . . . , (28) 
where = arcsin£. Thus, according to (|26[) the following dispersion relation is obtained for 7 = 5/3: 

A = A (m,fe) = m ± m = 0,±l,±2,.... (29) 

The admissible values of en, and are given by either (|27|) or (|28|l . These eigenvalues describe stable stratified sound waves 
that propagate independently on the surfaces of each individual cylindrical shell. In a co-rotating frame those are standing 
sound waves whose dependence on the axial coordinate is described by Eqs. (|27[) or (|28|) . where k plays the role of an axial 
wave number. Thus, according to (|27[) or (|28[) as k is increased the waves concentrate in growing numbers towards the disc's 
edges. It is finally emphasized that the compressibility effects in the vertical direction are indeed essential for the thin disc 
approximation as the sound vertical crossing time is of the order of fi -1 . 

The seemingly simple dispersion relations (|20[) and (|29[) actually entail an important conclusion. As they exhaust all 
possible eigenvalues of the linearized system it means that in the limit of thin disc approximation, the disc is spectrally stable 
against hydrodynamic perturbations. Thus, the asymptotic in time behavior of the linearized system exhibit two stable families 
of normal mode oscillations within the thin discs; one consists of vibrating rings due to inertia-coriolis effects, while the other 
one imposes density fluctuation on each individual cylindrical shell in the form of standing sound waves. In particular, it is 
obvious that no unstable strato-rotational modes may be excited in such systems. This is due to the absence of radial and 
azimuthal perturbed pressure gradients from Eqs. (|13[) and (|14[1 . which leads to the decoupling of the vertical density (strato) 
from the planar (rotational) dynamics. In contrast to that picture, in Couette flows, such pressure gradients play an important 
role in coupling of the planar and vertical perturbations and thus giving rise to the strato-rotational instability. The reason 
for the importance of the planar pressure gradients in Couette flows is the presence of a rigid wall, which is absent in rotating 
discs in thin disc approximation. 

Having clarified the time asymptotic (marginal) stability of rotating discs, this is by no means the end of the hydro- 
dynamical story. Short time dynamics may exhibit significant amplification of initial small perturbations. As will be shown 
below, two mechanisms are responsible for such behavior, namely resonant as well as non-resonant excitation of vertical sound 
waves by planar inertia-coriolis modes. Furthermore, the algebraically growing p ertur bations can create th e conditions for 



waves by planar mertia-conohs modes, furthermore, the algebraically growing perturbat ions can create th e conditions lor 
the re-appearing of the SRIs. It should be mentioned that Umurhan et al. ( 20061 ) and Rebusco et al. ( 20091 ) have described 



the non resonant mechanism in discs that are characterized by a phenomenological viscosity, by expansion in a single small 
parameter. The current work explores asymptotic expansions in the small but finite aspect ratio of the steady-state disc, e, for 
both equilibrium and perturbations, while the amplitude of the perturbations is assumed to be infinitesimal small independent 
value. The resonant mechanism is described here for the first time. 



4 NON-RESONANTLY DRIVEN SOUND WAVES 

Returning to Eqs. (|15|) and (|16[) that describe the dynamics of the vertical sound waves, it is clear that the latter may be 
excited by the inertia-coriolis modes that are eigen functions of the decoupled system (|13[) -()14 [ ). Moreover, due to the rotation 
shear, the amplitude of the driven sound modes grows linearly in time. This is illustrated by Eq. (|16[) which may be cast with 
the aid of expressions (I19|) - (|22|) in the following form of the auxiliary inhomogeneous equation: 



1 dn' dn dv' z ( ... dfi 1 dF 1 drfl\ , .m , . . 

n^T + 90 + 97 = ~r A 17+ f dF + rn^) Vr+l T v *- (30) 
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The solutions of the inhomogeneous Eqs. (|15[) - (|16l ) are represented now in the following way: 

{v' r , v'^, v z , n'} — exp(iut — im(f>){v r , v^, v z , h} + c.c. (31) 

Here cj = Afi(r) (A = const) represents the eigen-frequency of the inertia-coriolis waves, i.e., A = m ± 1, (m = 0, ±1, ±2, . . .) 
according to Eq. (|20|l . In addition it is assumed that 

v r = vi 0) , v$ = , v„ = v[ 0) + itnvi 1} , h = n (0) + itQ.h m . (32) 

Here the superscript denotes the order in time; Vr , are given by Eqs. (|2ip . (|22p while relations for n*-*' and vi (k = 0, 1) 
are derived after some algebra by substituting (|3ip and (|32|l into (|15[1 - (|16p and equating equal powers of t: 



±m n + __ = _ x __ Vr> (33) 



±^-^fn« + a| (C )^ = 0, (34) 
and 

± iftW + = -iftW - * + (35) 

±i ^ 0) - ^ti|^ (0) +c " (c) ^r = (36) 

with the following boundary conditions: 

n (0) = n (1) = for £ = ±1. (37) 

As is indeed seen from Eqs. (|33[) - ((36} the sound waves are driven by the inertial mode, while the linear growth in time 
is entirely due to effect of the rotation shear. Substituting now Eqs. pip . (|22p into Eqs. ((33} , ([34} and ((35} , ((36} yields the 
following two decoupled subsystems: 

<■ - ^ - + - ^m^- 

0<« = ±.(2 - 7K „<" ± ,1^1(1 - C 2 )^. (39) 
and 

,2,d 2 n^ 2 7 -3 dn<°> 3-7^0 4 „ m . 1 r F dCl . F 2 drFQi 

0<°> = ±i(2 - 7 K- (0) ±^(1 - C 2 )^ T (41) 

Equations (|38p and (|40p , subject to the boundary conditions (|37p and either symmetric or anti-symmetric scaling function 
G(£) are consistent with either symmetric or anti-symmetric solutions n^ k \Q (k = 0, 1), i.e: 

ftW(C) = ftW(-C), (42) 

ftW(C) = -ii W (-0- (43) 

A particular family of solutions is considered which is characterized by the following form of the free functions F(r) and 
G(C) for 7 = 5/3 (Appendix A): 

F(r)=r, G(C) = G 3 cos 3 9 + G 5 cos 5 6> + G 7 cos 7 0, (<9 = arcsinC). (44) 

As is shown in Appendix A3, such choice of F(r) and G(£) results in radial independence of the perturbed number density 
and its regularity at the disc edges. In that case the problem is reduced to two second order ordinary differential equations 
for n*- ' and n' 1 ' complemented by the corresponding boundary conditions ((42} or (|43[) at the disc edges, while ul ' and vi 1 ' 
may be calculated from (|39p and (|4ip . Thus, using (|44p the following symmetric solutions for rr ' and n' 1 ' are obtained: 

h {k) (0 = iVf 5 cos 3 6 + ivf 5 cos 5 6 + N$ k) cos 7 6, (45) 

where the coefficients AT^ and Gj (j — 3,5,7, k — 0, 1) are independent of the radius. Thus, for each value m there are two 
linearly growing modes except for m = ±1, where for each m corresponds only a single growing mode. It is finally noted that 
by the special construction of the solution in Appendix A3 both right hand sides of Eqs. (|38p and ((40} satisfy the solvability 
condition, namely that they are orthogonal to the solution of the common homogeneous part. 
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The special choice of the free function G(() presented in Eq. (|44f) indeed assures the regularity of all components of the 
velocity vector V'. However, such regularity requirements are rather too stringent and actually the sharp boundary model 
used in the current work does allow certain degree of singularity in the velocity. Thus, more relaxed conditions that are easily 
satisfied are that the perturbed components of the mass flux (or the momentum per unit volume), i.e. nV' is zero at the 
disc's edges. Indeed, such singularity in the velocity components is a direct consequence of the sharp boundary disc with 
zero equilibrium-density at the horizontal disc edges (£ = ±1) that has been employed in the present investigation. In fact, 
astrophysical discs have no well defined boundaries, the disc boundary is not sharp, and the vacuum boundary conditions 
imposed at a finite height are rather a good approximation than an exact model. According to a more realistic model that 



assumes a diffused disc with exponentially decreasing density [see e.g. iBalbus fc Hawlevl (|1991f ) for thin isothermal discs], the 
density has a small but finite value at the effective boundary of the disc. With that in mind the singularity at the disc edges 
may be regularized within a model with a small but non zero pressure value at the effective disc edges. The regularized model 
yields large but finite values of the fluid velocities at the disc edge that demonstrate that non modal growth may indeed leads 
to significant amplification of the perturbations. 
An example of a singular case is: 

F(r) = r, G(C) = Gi cos 9, (9 = arcsinC). (46) 

Again such choice of F(r) provides the solutions for (k = 0, 1) which are independent of the radius: 

n {k) = N^cosO, (N^ = T«|(m±l)Gi, N[ 0) = -i(4m±5)Gi). (47) 

In this case the velocity indeed has singularity at the disc's edges, however due to the more rapid decrease of the steady state 
number density the mass flux is zero. 



RESONANTLY DRIVEN SOUND WAVES 



As can be seen from Eq. (|30[) . stable sound modes may be excited resonantly by the inertial modes if any pair of respective 
eigen-values {20]) (with x = fi(r), uj = A^ m) fi(r), where A^ m) = mil) and ([29]) are equal, i.e. A^ m) = \ ( ±- k) (m = 0, ±1, ±2, . . ., 
k — 0, 1). It is easy to show that this happens only for the first anti-symmetric mode (|28[) for arbitrary m (7 = 5/3) and: 

k = ±l, a k = ±2, A^ m ' fe) =m±l, h = ±N 2 sin 29, (9 = arcsin <) . (48) 

The solution of Eq. (|30|) then, that describes the resonant coupling between sound and inertial waves in Keplerian discs 
is thus represented by: 

{v' r , w^, v' z , n'} = exp(iut — im(j)){v r , v^, v z , n} + ex., oj — XQ(r), (A = const). (49) 
The combination of resonant coupling with the effect of the rotation shear leads to the following form of the solution: 

v r = e\ = vf\ V* = vi 0) + itSivP + l!M!fi( 2 >, n = n<°> + itQn^ + ^!ft( 2 ) , (50) 

where the superscripts denote the order in time; i>l°\ v^' are given by (|21|) . (|22l) while vi k \ nS k \ (k = 0, 1, 2) are determined 
by the relations below. It should be noted that resonant interaction between the two modes of wave propagation may occur 
also in a shearless rotation. In that case the amplitude of the sound waves grows only linearly in time. It is the combined 
effect of resonant interaction and rotation shear that gives rise to a quadratic growth in time that is described in Eq. (I50|) . 
Substituting (|49|> - (|50[) into (1131) ) - (|16[) and equating equal powers oft yields 

± «W + ^ = , (51) 

±«. W "^f * W +£CO^=0, (52) 
and 

dvP _( 2 ) . 1 dCl „ . . 

i m H = -iw - A — — v r (53) 

oQ il ar 

± tv W ^-7T-n m +c% (O-apr- =0, 54 

7 — 1 aQ oQ 

and 

± ift (0) + MP = -±*™L+ ™v„ (55) 

d(, r\l or r 
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± ^ - ^f^ 0) + ^ (0^ - (56) 

while the planar dynamics is described by Eqs. (|21|l . ()22[) for v r and v^. The number densities (k = 0, 1, 2) are supplemented 
by the boundary conditions: 

h m = n (1) = n (2) = for C = ±1- (57) 
Equations (1551) and (|56|) for n^ 2 ' and form a homogeneous sub-system for the sound mode: 

0?' = ±.(2 - 7 )C»™ ± <3L=i(l - C ")^. ( 59 , 

The only solutions of Eq. ()58[) that satisfy the resonance conditions and which are zero at £ = ±1 and satisfy (|57|) are given 
by ( 7 = 5/3): 

n (2) = ±iV 2 (2) sin 20, vf ] = iAT 2 (2) [(2 - 7) sin 9 sin 20 + (7 - 1) cos cos 20], (0 = arcsinC)- (60) 

In general, A^ may depend on the radial variable, however, as an example it will be considered as a constant. Also, it is 
noted that in contrast to the non-resonant case, here singularity in the velocity is inevitable. However, as stated above, the 
less restrictive condition is satisfied as the mass flux is regular and vanishes at the disc's edges. 
Moving on the following equations for n' 1 ', wi 1 ' and n' ', v^: 

(1 - C 2 )^ - ^2 C ^ + 2^n« = T ^-n< 2 > =p 3^6(0, (61) 
aC, 2 7 — 1 aQ 7— 1 7—1 7 — 1 

= ±i(2 - 7 )C- (1) ± " C 2 )^, (62) 

and 

(1 - < 2 )^ - ^2 C ^ + 2^n<°) = -^n< 2 > =p -^n«" - 2i^±G, (63) 
a(, 2 7 — 1 a(, 7 — 1 7 — 1 7 — 1 7 — 1 

4 0) = ±i(2 - 7 )C- (0) ± - C 2 )^ T 0«. (64) 

The homogeneous part of Eq. (|6ip complemented by corresponding boundary conditions for rr 1 ' coincides with the above 
eigenvalue problem for fi' 2 '. Hence non-trivial solutions of the inhomogeneous equations exist if the right hand sides of (|61[) 
satisfies the solvability condition, namely, that it is orthogonal to hi as it is given in (|60|) . This means that in order to obtain 
resonant coupling, the vertical dependence of the inertia-coriolis waves must contain a term of the following shape: 

G(C) = G 2 sin 20, (0 = arcsin <) , (65) 



where G2 is an arbitrary amplitude. For simplicity it will be assumed that Eq. (|65[) indeed describes the entire profile of the 

driving waves and in addition that F(r) = r. As will be seen later, the latter results in simple radial independent solutions. 

(2) 

Imposing now the solvability condition, A^ is given by: 

=Ti|(n»±l)G a . (66) 

Eq. (|61[) for nS 1 ' may be solved now. Due to imposing the solvability condition and the special form of G given in Eq. 
(|65[) . the right hand side of that equation is zero. Hence its solution is given by A^ 1 ' sin 20. Turning now to Eq. (|63[) it is 
immediately seen that the same solvability condition determines now AT 2 (1) to be: 

^2 (1) =±ij(m=Fl)G 2 . (67) 



Summarizing, the solutions for n*- 1 ', 0* and n/ ', vi * 1 are (0 = arcsin Q: 



h w = AT 2 (1) sin20, = ±iAT 2 C1) [(2 - 7) sin sin 20 + (7-1) cos cos 20], (68) 

n (0) = A"^ O) sin20, v ( z 0) = i[±N^ 0) - jV 2 (1) n(r)][(2 - 7) sin sin 20 + (7 - 1) cos cos 20]. (69) 

In the above solution (|66[) - (|69[) there are two arbitrary constants G2 and A^ ' . The first one represents the arbitrary amplitude 
of the driving inertia-coriolis waves while the second one describes free sound oscillations. Finally, it is noticed again that 
while the vertical component of the perturbed velocity is singular at the disc's edges, the corresponding component of the 
mass flux tends to zero. 
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6 CONCLUSIONS 



A comprehensive study of the long as well as short time hydrodynamic response of Keplerian rotating discs to small pertur- 
bations has been carried out for realistic thin disc geometry. Explicit analytical solutions of the classical eigenvalue problem 
reveal that thin discs are hydrodynamically marginally stable asymptotically in time. Such a picture naturally leads to the 
ruling out of the strato-rotational instabilities as primary contributors to the dynamical evolution of thin ultrasonically rotat- 
ing discs. This has been shown to be a direct result of the decoupling of the planar dynamics from the vertical acoustics, which 
in turn is a result of the negligible role played by radial pressure gradients. The latter is indeed the agent of coupling of the 
strato and rotational modes that lead to instability in the case of Couette flows due to the presence of solid radial boundaries. 
Notwithstanding the spectral stability of thin rotating discs, they have been shown to exhibit intense hydrodynamical activity 
in the form of algebraic short-term growth of initially small perturbations. Two mechanisms have been identified and studied 
that lead to such non modal dynamical behavior. The first one is the non resonant excitation of co-rotating standing sound 
waves on individual cylindrical shells by planar inertia-coriolis modes. Such coupling between the modes occurs due to the 
rotation shear and results in a linear in time growth of the perturbations. The second mechanism is the resonant driving of 
the sound waves by the inertia-coriolis modes. While non resonant coupling relies on the non normal nature of the linearized 
system of the dynamical equations, resonant coupling exists also for a shearless rotation. However, in the presence of rota- 
tion shear the perturbations grow quadratically in time. Both mechanisms are powered by the compressible motion of the 
fluid. Thus, compressible inertia-coriolis stable oscillations continuously pump energy from the sheared steady state flow and 
transfer it (resonantly as well as non resonantly) to the continuously algebraically growing sound waves. Compressibility is 
indeed inevitable due to the supersonic rotation of the steady-state disc combined with its small vertical dimensions which 
make the vertical sound crossing time of the order of a rotation period. Simple explicit analytical solutions were obtained for 
all possible cases. It should be stressed that non modal growth has been exhibited for all admissible parameters of the system. 

Along side with the non modal growth of small perturbations, the rotation shear is also responsible for the linear growth 
in time of an effective radial wave number. Such a process naturally gives rise to enhanced perturbed radial pressure gradients. 
In light of the discussion of the special role played by the latter in the onset of the SRI, it may be speculated, that such 
a development may reinstate the SRI as an important dynamical process due to secondary instabilities that occur in the 
wake of the primary non modal hydrodynamical growth. Thus, the linear as well as quadratic amplification of initial small 
perturbation as well as the secondary instabilities they entail could provide a viable mechanism for sustaining turbulence and 
consequently to enhanced transport coefficients. 

A detailed description of nonmodal growth has been presented before in lUmurhan et al. | |200ct and lRebusco et alj (|2009h ■ 
The main differences between those works and the current investigation are: 

• In Umurhan et al. ( 20061 ) and Rebusco et al. (2009) a finite amplitude non-linear solution to the problem of the dynamical 
evolution of the full equations was sought. In that approach it is assumed that solutions of the equations, both steady as well 
as dynamical, may be expanded in a single small parameter, namely, e. Solutions are then determined by iteratively solving 
the resulting equations at each order of e. Consequently, the forcing terms at each order are the solutions of the previous order 
which results in algebraic growth. In contrast, in the current work a classical linearization procedure is employed in which 
an extra independent small parameter is introduced that measures the amplitude of the initial perturbations. It is due to 
the negligible role played by the radial pressure gradients that the resulting linear system of equations is decoupled into two 
subsystems, one of which describes the perturbed planar velocities while the other one determines the dynamical evolution 
of the vertical sound waves. Moreover, the planar modes act as forcing terms for the vertical sound waves. This mechanism 
give s rise to the resonant a s wel l as the nonresonant a lgebraic growth of the latter. 

• llJmurhan et al] <|2006h and iRebusco et"aH (|2009h have demonstrated the occurrence of linear growth (due to the shear 



of the steady state rotation) that results from the nonresonant driving of the vertical sound waves by the planar dynamics. 
In the current work, in addition to that linear growth, a resonant process is identified which results in a quadratic growth in 
time. Such second order growth is a result of the combined effect of the resonance between the vertical and and the planar 
modes and the shear of the steady state rotation. Furthermore, the resonance is a global one in the sense that it occurs all 

over the disk; the planar vibrations of a ring of any radius resonantly ratt le the cylindrica l shell where it resides. 

• Unlike the inviscid limit that is investigated in the current work, lUmurhan et al. (|2006h and IRebusco et alj (2009) 
present the solution o f the viscous problem whe re the viscous stresses are modelled by the Shakura-Sunyaev a parameter. 
Umurhan et al. ( 2006h and Rebusco et al. (2009) refer to values around a — 10 3 as representing "realis t ic" va lues as they 
are close to the result found numerically in subcritical hydrodynamic transitions in iLesur fc Longaretti (2005). Numerical 
solutions indicate that for such small values of the viscosity (or large values of the Reynold number), the algebraic growth 
proceeds undisturbed at least for 10 3 rotation times, during which the amplitude grows by a factor of I0 4 , before viscous 
damping starts to play any significant role. Thus, the dynamical evolution during the first non-viscous long period of time is 
indepe ndent of the viscosity and is well described by the inviscid limit. Indeed, the eigenvalues of the planar dynam i cs [Eq . 
(63) in Rebusco et al. ( 20091 )] as well as the temporal behavior of the driven sound waves [Eq. (66) in Rebusco et al. ( 20091 )] 
transit to Eq. (|20[) and Eqs. (|3f | Eqs. ((32} , respectively of the current work in the limit a — > 0. The singular limit of zero 
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viscosity is revealed then only in the spatial profiles of the eigenfunctions. Thus, the spatial profiles obtained in the current 
work may be considered as outer solutions of a boundary layer problem which are valid in most parts of the disk except 
perhaps within a thin region next to the vertical edges of the disk. The advantage of the inviscid approach however is that it 
enables to derive simple analytical expressions for the eigenvalues and the eigenfunctions such as those giv en in Eqs. (1271) - (1291) . 
and co nsequently to identify the resonant quadratic growth. The inviscid problem has also been solved bv lUmurhan fc Shaviv 
( 2005h where they have described the nonresonant growth with a single parameter expansion (see first item). 



Obviously, the algebraically growing perturbations (both resonant as well as non resonant) do not keep growing without 
bound. Thus, nonlinear effects as well as v iscous damping are e xpected to quench the linear or quadratic growth of the 
perturbations. As has been demonstrated in iRebusco etail (|2009l ). for realistic values of the viscosity, the algebraic growth 
of nonresonantly driven sound waves takes place unhindered for at least 10 3 rotation times before viscous damping kicks in. 
During that time the initially small perturbations grow significantly by few orders of magnitude. It is expected therefore, that 
nonlinear redistribution of the perturbations energy among other length scales will occur before the perturbations are damped 
by viscous effects. In the case of resonantly driven sound waves the growth of the perturbation is even faster. In addition, as 
\k r \ ~ \t\dVl/dr\ (see Section 3.1) the time needed to develop significant perturbed radial pressure gradients is of the order of 
(Ae) _1 . According to lRebusco et al. I |2009h there is enough time for that to occur before viscous effects become important and 
even before nonlinear effects take over. Moreo ver, the time to e stabli sh si gnificant radial gradien t s decr eases as the azimuthal 
mode number is increased. Also, according to lDubrulle et al. I (|2005h and IShalvbkov fc Ritdigerl (|2005l ). the maximal growth 
rates of the SRI are achieved for Froude numbers of the order of ~ 1 — 10 and are of the order of 10 _1 inverse rotation 
periods. This, again, leaves ample time for the growth of secondary SRI i nstab ilities before t he perturbations are k i lled o ff by 
viscous stresses. It should be mentioned however that in Dubrulle et al. (2005) as well as in lShalvbkov fc Ritdigerl (|2005T ) (as 
indeed in most SRI investigations) the Froude number has been assumed to be constant. This is very different from realistic 
thin disk configurations where the Froude number is given according to eqs. ([Toll and (H2} by Fr 2 = 3C 2 /(1 - C 2 ) and thus 
varies between zero and infinity over the small vertical extent of the disk. Hence, the numbers cited above are used just as a 
crude estimation and the true nature of a possible secondary instability may be investigated at this stage only by numerical 
simulations. 

It should finally be stressed once again that the calculations presented above are valid in the Keplerian portion of the disk, 
namely far away from the radius of zero torque. In regions around the latter, which are much closer to the central object, the 
radial pressure gradients may play an important role in the dynamics of small perturbations, as is shown in lRevnolds fc Miller 
[ 2009h who studied numerically g-mode trapping that occur within 8-32 Schwarzschild radii from a central black hole. 



ACKNOWLEDGMENTS 

This work has been carried out while one of us (M.M) was a Mercator Professor at the university of Potsdam. He would 
like to thank the Physics department at the University of Potsdam and the Astrophysikalisches Institut Potsdam for their 
hospitality. The authors thank Georgy Burde for useful and insightful discussions. 



APPENDIX A: 

Al The basic boundary-value problem 

Equations (|26p . (|38p . (|40|l . (|58[1 and (|61[l . (|63|l have the following general form (it was chosen F(r) = r everywhere through 
that section): 

(! ~ ^ ~ + ^ = C K X ' ° 2 = ^T [(A - -) 2 + 2 - 7l- (Al) 

Equation (|A1|) is supplemented by the boundary conditions at the disc edges. In addition, the solutions of Eq. (|A1[) satisfy the 
same symmetry as that of the driving term on the right hand side, namely either symmetric (^(C) = C)) or anti-symmetry 
(*(C) = -*(— 0) conditions: 
Hn 

— = for C = 0, n = for ( = 1, (A2) 
aQ 

n = for C = 0, n = for C = 1- (A3) 
Transforming to a new independent variable 6 = arcsin£, problem (|A1[) - (|A3[) assumes the following form: 

__-2__tan0— + a 2 n = *(C), -rr/2 < G < tt/2, (A4) 
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fjr) 

— = for 9 = 0, n = for 9 = vr/2, (A5) 
n = for 9 = 0, n = for = vr/2. (A6) 



For arbitrary values of 7 the solutions of Eq. (|A4[) may be expressed in terms of the hypergeometric functions. However, for 
the specific value 7 = 5/3 Eq. (1A4[) becomes remarkably simple. Conveniently therefore, considering from now on the case 
with 7 = 5/3 as well as Keplerian rotation, = r~ 3 ^ 2 . 

A2 The eigenvalue problem for the sound mode 

The eigenvalue problem for the sound modes is given by (see Eq. (|26[) for 7 = 5/3): 

Jl_ +Q 2 n = 0, O<0<tt/2, (A7) 

n = for 9 = ±tt/2. (A8) 
The general solution of Eq. (|A7|I is: 

n = Ci cos a9 + C 2 sin a0, ( A9) 



where cos a9 and sin a9 are the symmetric and anti-symmetric eigenfunctions. Imposing boundary conditions (|A5[l and (|A6[> . 
the following two families of eigenvalues and eigenfunctions obtained: 

a = 2k + l, n = N 2k +i cos(2fc + 1)0, fc = 0, ±1, ±2, . . . , (A10) 
a = 2k, n = N 2k sin 2fc0, fc = ±l,±2, .... (All) 

A3 The regular symmetric solution for the non-resonantly driven sound mode 

The following choice of a symmetric driving function *&(C) in Eq. (|A4|I is considered (see Eqs. (|38[1 and (|40[) with 7 = 5/3): 

]^+n=^(0 = ^ 3 co S s 9 + ^ 5 cos 5 9 + ^ 7 cos 7 9, < 9 < tt/2, (A12) 
4 a0 2 4 

— = for 9 = 0, n = for = tt/2, (A13) 
a0 

where A = m ± 1, a 2 = 4. The motivation for such a choice of the right hand side in eq. ()A12fl is the following: recalling 
solutions (|20[) and (|21|) for and vj,, as well as definition (|19[) . and the steady state solution for the number density (|12|l 
[n ~ (1 — ^ 2 ) 3 / 2 ~ cos 3 9 ], it is seen that the solution for V r ' and that are regular at the disk edges £ 2 = 1 are given by 
a function G that contains powers of third order and higher of cos#. It will consequently be shown in this subsection that 
similar regularity requirements on VL imply that G must contain at least also the next two higher odd powers of cos 9 whose 
amplitudes are determined by the arbitrary amplitude of the cubic term. 
The general solution of (IA12I) ') for symmetric modes is 

n = C 2 cos 29 + Ci cos 9 + C 3 cos 30 + C 5 cos 50 + C 7 cos 79, (A14) 

where from (|A13[) at 9 — tt/2 follows that C 2 = 0, and 

Ci=Vs + |v6 + §^, C 3 = -~V>3 - ~V>5 - §V7, C 5 = ~i*-±i>T, C 7 = -±1». (A15) 
For further convenience (|A14p is rewritten in the form 

n = Ni cos 9 + N 3 cos 3 + N 5 cos 5 + N 7 cos 7 0, (A16) 
where 

Ni = Ci-3C 3 +5C 5 -7C 7 , Af 3 = 4C3-20C5+56C7, A^s = I6C5-H2C7, N 7 = 64C 7 , 
or using (A16) 

„ 8, 32, 64, lr 4 , 16 , 32 , 4 , 8 , 4 , , A1 _ 

Ni = -ipz + —ip 5 + —ip 7 , N 3 = --ip 3 -—ij 5 -—ip 7 , N s = -— V> 7 , A 7 = -— Vr- (A17) 



Following the regularity requirement on and Eqs. (|39[) the coefficient of cos9 ~ (1 — (J 2 ) 1 / 2 in (|A16() must is zero, namely: 

8 , 32 , 64 , 
= 5^+21^5 + ^ = 0. 

Applying relations (|A12|) - (|A18|I to Eq. J3S]| for yields: 



iVi = |Va + + ^7 = 0. (A18) 
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G = G 3 cos 3 6 + G 5 cos 5 6 + G 7 cos 7 0, (A19) 
and 



= T^AG.,, A = m ± 1. (A20) 
8 

Here ^/>j are the coefficient for the right hand side of Eq. (|38l) for ft' 1 ', the coefficients G-, will be specified below (j = 3, 5, 7). 
The above relations determine ft' 1 ' 

ft (1) = n[ 1] cos 6 + iVg 1 ' cos 3 6> + A^ 1 ' cos 5 6 + cos 7 6, (A21) 
where JV- 1 ' are expressed in terms of ipj 1 ^, (j = 1, 3, 5, 7), through relations ()A17|) 

NP =TiA(^G 3 + yG 5 + ^G 7 ), iVf =±iA(^G 3 + ^G5 + ^G 7 ), A^ 5 (1) = ±iA(^G 5 + ~G 7 ), A< x) = ±iA^G 7 , (A22) 
with the regularity condition 

Af' = =rA(|g 3 + -^G 5 + f G 7 ) = 0, A = m ± 1. (A23) 
5 7 5 

Hence ft' 1 ' is determined through Gj (j = 3, 5, 7), which will be specified below. 
Applying relations (|A12|I - (|A18[) to Eq. (gDJ) for n (0) yields: 
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^Pf =t\n^ -\i{m±2)G 3 , (j = 1,3,5,7). (A24) 



Substituting (|A22|I into (]A24|) yields 

^(°) = --i( m ±2)Gi, ^3°' = -— i(Um±19)G 3 ~ i(m±l)G 6 --i(m±l)G r , 
4 20 T 5 

4 0) = -^i(10m±17)G 5 -^i(m±l)G7, 4 0) = -^<(6m± H)G r , (A25) 

and 

ft (0) = ATf cos 61 + AT 3 (0) cos 3 6> + A^ 0) cos 5 6 + cos 7 9. (A26) 

The coefficients A?"j ' are expressed through ijjj according to (|A17[I and then by using (IA251 through Gj, (j = 1,3,5,7). 
Once again, requiring the regularity of Vl and employing Eq. (|41[) results in: 

N[ 0) = -i— (14m±19)G 3 - i— (113m ± 148)G 5 - i — (40m ± 37)G 7 = 0. (A27) 
25 245 175 

Finally, setting, e.g. G3 = 1 as an arbitrary amplitude of the inertia-coriolis driving waves, Eqs. (|A23|) and (IA27[) determine 

the coefficients G5 and G 7 . Since both F(r) and G(£) are known, the corresponding velocities v r , V4,, Vi 1 [k = 0, 1) may be 

calculated along with the number density. 
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